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ABSTRACT 

We report on observations of the unusual neutron-star binary system 
FIRST J102347.6+003841 carried out using the XMM-Newton satelhte. This system con- 
sists of a radio milhsecond pulsar in an 0.198-day orbit with a ~ 0.2 Mq Roche-lobe- filling 
companion, and appears to have had an accretion disk in 2001. We observe a hard power-law 
spectrum (F — 1.26(4)) with a possible thermal component, and orbital variability in X-ray flux 
and possibly hardness of the X-rays. We also detect probable pulsations at the pulsar period 
(single-trial significance ~ 4.5<t from an 11(2)% modulation), which would make this the first 
system in which both orbital and rotational X-ray pulsations arc detected. We interpret the 
emission as a combination of X-rays from the pulsar itself and from a shock where material 
overflowing the companion meets the pulsar wind. The similarity of this X-ray emission to that 
seen from other millisecond pulsar binary systems, in particular 47 Tuc W (PSR J0024— 7204W) 
and PSR J1740— 5340, suggests that they may also undergo disk episodes similar to that seen in 
J1023 in 2001. 

Subject headings: pulsars: general — pulsars: individual (PSR J1023-I-0038) — stars: neutron — X-rays: 
stars 
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1. Introduction 

Radio millisecond pulsars are thought to have 
been spun up to their current rapid rotation by a 
long period of accretion in a Low-Mass X-ray Bi- 
nary (LMXB) . This general idea is well supported 
by many lines of evidence, but the transition be- 
tween LMXB and radio millisecond pulsar remains 
mysterious. In particular, models of binary evolu- 
tion often predict that mass transfer should con- 



tinue indefinitely (Deloye 2008, ), but ionized ma- 
terial entering a neutron star's magnetosphere is 
expected to "short out" magnetospheric emission 



( Stella et al.||1994^ . How then can we explain the 
observed radio millisecond pulsars? 

An important system for studying this question 
was recently discovered: FIRST J102347.6-h003841 
(hereafter J1023) is a 0.198-day binary containing 
a 1.69-ms radio millisecond pulsar, PSR J1023-I-0038. 
This system is unique in that it shows evidence for 
having had an accretion disk as recently as 2001, 



but is now active as a radio pulsar (Archibald 



et al. 2009 Wang et al. 2009). This suggests 



that we are observing the end of an LMXB phase 
and the birth of a radio millisecond pulsar, dur- 
ing which the system undergoes transient LMXB 
phases. This is consistent with a scenario put for- 



ward by Burderi et al. (2001) in which LMXBs, 



late in their evolution, may have episodes of mass 
transfer but spend most of their time in "radio 
ejection" phases, during which the wind from an 
active radio MSP sweeps material overflowing the 
companion out of the system. 

In its current quiescent state, J1023 has been 
extensively observed with optical and radio tele- 
scopes. The millisecond pulsar's companion ap- 
pears as a F ^ 17 magnitude star with a G- 
like absorption-line spectrum. The brightness and 
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temperature show a 0.198-day modulation consis- 
tent with heating due to an irradiation luminosity 
of roughly 2Lq from the primary (Thorstensen 



& Armstrong 2005). The combination of pulsar 



timing and optical radial velocity measurements 
shows that the mass ratio is 7.1(1), while assuming 
the pulsar mass is ~ 1.4M0 implies an inclination 
angle of ~ 46°, a companion mass of ~ 0.2 Mq 
and, based on Roche lobe size and optical appar- 
ent magnitude, a distance of '-^ 1.3 kpc (Archibald 



et al. 20091. At radio wavelengths, the system 



shows variable, frequency-dependent eclipses near 
the orbital phase at which the companion is closest 
to our line of sight to the pulsar, though the com- 
panion's Roche lobe never actually intersects this 
line of sight. Irregular short eclipses are also ob- 
served at all phases ( Archibald et al.|["2009 |. Such 
variable radio eclipses have been observed in other 
millisecond pulsar systems, both in the so-called 
"black widow" systems in which the companion 
is thought to be a very low-mass degenerate star 



being ablated by the pulsar wind ( Fruchter et al. 
1988 Stappers et al.|1996 ), and in several systems 



more similar to J1023, with companions that ap- 



pear to be unevolved and Roche-lobe-filling ( Freire 
et al.||2003[ [D'Amico et al.||2001[ ). Comparison of 



J1023 with these systems, as well as with LMXBs 
in quiescence, could clarify how accretion comes 
to an end and radio pulsars become active. 

Prior to the discovery that J1023 contained 
a millisecond pulsar, and unfortunately after it 
had returned to quiescence, 'Homer et all (|2006|) 



observed J1023 with XMM-Newton. They con- 
cluded that it had a hard power-law spectrum, 
variability not necessarily correlated with the or- 
bital phase, and a 0.5-10 keV X-ray luminosity 
of 2.5 X 10^2 gi-g g-i (O.O65L0) at a distance, 
based on estimates available at the time, of 2 kpc. 
This supported the suggestion of [Thorstensen &:] 



Armstrong (2005) that J1023 had a NS primary 



rather than a white dwarf as had previously been 
thought. 

Not only does the detection of J1023 as a radio 
millisecond pulsar confirm that it is a NS binary, 
radio timing provides an ephemeris predicting the 
rotation of the pulsar. For this reason, we ob- 
served J1023 for a longer time and with high time 
resolution using XMM-Newton. Our observations 
confirm the spectral observations of [Homer et al. 
(2006), show evidence for orbital variability, and 
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show for the first time probable X-ray pulsations 
at the pulsar period. 

In section 2 we summarize the observations we 
used and our initial data reduction. In section 
3 we describe and give results from our spectral 
analysis and our search for orbital and rotational 
modulation in the X-ray data. In section 4 we 
discuss the implications for the nature of J1023's 
X-ray emission and the relation of J1023 to sim- 
ilar systems, and in section 5 we summarize and 
suggest further avenues of researcli. 

2. Observations 

The European Photon Imaging Camera, EPIC, 
aboard the XMM-Newton satellite, has two MOS 
cameras ( Turner et al.|200r ) and one PN ( Striider 
et al.|2001 ). All three cameras may be operated in 
an imaging mode, but the PN camera can also be 
run in a fast timing mode, in which the photons 
from one CCD are accumulated along one dimen- 
sion, losing two-dimensional imaging, but achiev- 
ing a time resolution of 29 fj,s. Calibration uncer- 
tainties limit the absolute time accuracy of XMM- 



cal modulation consistent with that reported by 



Thorstensen & Armstrong (2005), so we did no 



Newton to approximately 100 ^s ( Guainazzi et al 



2010). XMM-Newton also has a coaxial optical 



telescope, the Optical Monitor, and two reflection 
grating spectrometers, however J1023 is too faint 
to permit interesting studies with the latter. 

We obtained an XMM-Newton observation 
(OBSID 0560180801) of J1023 of 33 ks, corre- 
sponding to 1.94 binary orbits, on 2008 Nov 26. 
We operated the PN camera in fast timing mode 
and the MOS cameras in "prime full window" 
(imaging) mode, and we used the thin filter on all 
the EPIC cameras. The data were free of soft pro- 
ton flares, so we were able to use the entire expo- 
sure time. We also retrieved the data set (OBSID 
0203050201) used by [Homer et"al] ( [2006| from the 
XMM-Newton archive; it consists of a single 16 ks 
exposure (0.92 binary orbits) obtained on 2004 
May 12, with all EPIC instruments using "prime 
full window" mode and the thin filter. This data 
set is also free of soft proton flares, but since all 
three EPIC instruments were operated in imaging 
mode, in this case high-resolution timing data are 
not available. 

The Optical Monitor was operated in both ob- 
servations with the B filter and a mode with 10-s 
time resolution. Pipeline processing showed opti- 



further analysis of the optical data. 

We analyzed the X-ray data with the XMM- 
Newton Science Analysis Softwar^ version 9.0.C[^ 
and xspec version 12.5.0ac[^ 

We reprocessed the MOS and PN data with the 
emchain and epchain pipelines, respectively, then 
extracted photons meeting the recommended pat- 
tern, pulse invariant, and flag criteria. In particu- 
lar, this restricted MOS photons to the 0.2-12 keV 
energy range, PN imaging photons to the 0.13- 
15 keV energy range, and PN timing photons to 
the 0.25-15 keV energy range. For all analyses ex- 
cept the spectral fitting, we restricted the imaging 
photons from either type of camera to 0.2-10 keV 
for consistency, since this involved discarding only 
a handful of photons. 

The default pipeline processing showed no ev- 
idence for extended emission from J1023, so we 
used the SAS tool eregionanalyse to select ex- 
traction regions to optimize the ratio of source 
photons to the square root of source plus back- 
ground photons. These regions were circles of ra- 
dius 35" (enclosing about 90% and 85% of the 
1.5 keV flux on the MOS and PN cameras, respec- 
tively) in imaging observations. For the timing- 
mode PN observation we adjusted the extraction 
region by hand to optimize detection significance 
(see below), selecting a band 12" wide (enclos- 
ing about 60% of the 1.5 keV fiux). We selected 
background regions from nearby large circles in 
the imaging-mode observations, and from a band 
50" wide, to one side of the source, in the timing- 
mode PN data. All together, this yielded approx- 
imately 5700 background-subtracted source pho- 
tons in imaging mode and approximately 2400 
background-subtracted source photons in timing 
mode. 

We computed photon arrival times using the 
SAS barycentering tool, barycen, the precise op- 
tical position 10'^23™47.67^ -t-0°38^41.2" (J2000) 



given in the NOMAD catalog ( |Zacharias et al. 



^The XMM-Newton SAS is developed and maintained by 
the Science Operations Centre at the European Space As- 
tronomy Centre and the Survey Science Centre at the Uni- 
versity o f Leicester. 

%ttp : //xmm. esac ■ esa. lnt/sas/| 



http : //heasarc . nasa . gov/docs/xanadu/xspec/ 
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2004), and the DE405 solar system ephemeris. 



To produce pulsar rotational phases during the 
2008 observation, we used the tool tempcF] and 
the ephemeris given in Archibald et al. (p009). 



which is based on radio pulsar timing data brack- 
eting the X-ray observation epoch. The extremely 
sensitive phase-coherent timing procedure used by 
[Archibald et al.l revealed minuscule variations in 
the orbital period of J1023. These are modelled 
in the published ephemeris by expressing the or- 
bital period as a quadratic polynomial obtained 
by a fit to several months of 2008 timing data. 
Extrapolating this polynomial back to the 2004 
XMM-Newton observation results in an unreason- 
able and poorly constrained orbital period at the 
time of ^ —6 days, a clear indication that the 
orbital period variations are not actually given by 
such a quadratic polynomial. While it is necessary 
to model the orbital period derivatives to obtain 
good quality pulsar phase predictions, it is possi- 
ble to obtain adequate orbital phase predictions by 
using a simple model with constant orbital period 
of 0.19809620 days and (pulsar) ascending node of 
MJD 54801.970652993, based on the model from 



Archibald et al. ( 2009 ). This model gives adequate 
orbital phase predictions back to 2004, match- 



ing the orbital phases observed in Thorstensen 



& Armstrong (2005) to within 1% of a period as 



well as matching the ephemeris given in | Archibald] 
et al.l (|2009|) exactly on 2008 Dec 1. We there- 



fore used this model to evaluate the orbital phase 
of each barycentered photon in both our observa- 
tions. 

3. Results 

3.1. Spectral fitting 

We used xspec to fit each of several models to 
all EPIC imaging data sets from both epochs si- 
multaneously. To verify that the data sets were 
compatible, we also fit different absorbed power 
laws to the two observations, obtaining compati- 
ble spectral indices and normalizations. For the 
MOS cameras we restricted photons to the recom- 
mended energy range, namely 0.2-10 keV, and for 
the PN camera we used the recommended lower 
limit of 0.13 keV but used the upper limit of 
10 keV as there were not enough photons above 
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Fig. 1. — Absorbed neutron-star atmosphere plus 
power-law spectral fit to all EPIC imaging-mode 
photons from both observations. Each data set 
is plotted in a different color: black is the 2004 
MOSl, red the 2004 M0S2, green the 2004 PN, 
blue the 2008 MOSl, and cyan the 2008 M0S2 
data set. The vertical axis is normalized by the 
effective area, and the expected responses from 
the two additive components are overplotted sep- 
arately (the harder component is the power law). 
The lower panel is the ratio of model predicted flux 
to observed flux in each group of detector chan- 
nels. 



http : //www . atnf . csiro . au/research/pulsar/tempo/ 



4 



Table 1 

Spectral fits to all EPIC image data 





power-law 


neutron star atmosphere 
plus power law 


TVh (cni-2) 


< 5 X 10^9 


< 1 X 10^0 


photon index 


1.26(4) 


0.99(11) 


kT (keV) / photon index (F) 




0.12(2) 


Thermal emission radius (km)'' 






0.5-10 keV unabsorbed flux (erg cm^^ s^^) 


4.66(17) X 10-13 


4.9(3) X 10-" 


0.5-10 keV luminosity'' (erg s^^) 


9.4(4) X 10^1 


9.9(5) X 10^1 


Thermal fraction 




0.06(2) 


/degrees of freedom 


230/214 


208/212 


Null hypothesis probability 


0.21 


0.56 



Assuming a distance of 1.3 kpc. 



this to provide any constraint. We grouped the 
photons to obtain at least 20 counts per spectral 
bin so as to have approximately Gaussian statis- 
tics. 

We obtained an adequate fit to all imaging data 
sets with a power-law model including photoelec- 
tric absorption. We did not obtain adequate fits 
with an absorbed single black-body (x^ — 1203.7 
with 213 degrees of freedom, for a null hypoth- 
esis probability of ~ 10"^^^) or neutron-star at- 
mosphere model (x^ = 1007.1 with 213 degrees 
of freedom, for a null hypothesis probability of 

~ 10-103). 

The simple power-law model, with a photon in- 
dex of 1.26(4), already has a null hypothesis prob- 
ability of 0.21, which is satisfactory from a statis- 
tical point of view. Although the statistics do not 
call for a thermal component, we may expect one 
on physical grounds (for example from the surface 
of the neutron star) . We therefore also consider an 
absorbed neutron-star atmosphere plus power-law 
model to fit these data, as shown in Figure [l] This 
model, specified in xspec as phabs (nsa+pow) , is 
based on Zavlin et al. (1996). During fitting we 
froze the neutron star mass at I.4M0 and the ra- 
dius at 10 km for the purposes of redshift and light 
bending; the smaller effective thermal emission ra- 
dius suggests that the emission comes from a small 
"hot spot" on the surface. In light of the fact tha t 
the estimated S < 10** G ( [Archibald et al.||2009| ). 



we also selected a model in which the magnetic 
field has negligible effects on the atmosphere, i.e. 
< 10^ G. 

These two models are summarized in Table [T] 
Uncertainties given are 90% intervals returned 
by xspec's error command; luminosities are ob- 
tained by using the cf lux model component to es- 
timate unabsorbed and thermal fluxes in the 0.5- 
10 keV range and then using an assumed distance 
of 1.3 kpc to determine a luminosity. "Thermal 
fraction" is the fraction of this luminosity due to 
the thermal component of the spectrum, if any. 

In both cases, fitting for photoelectric absorp- 
tion in the model gave a very low upper bound 
on the neutral hydrogen column density, as noted 
in Homer et al. (2006). Such a low value is to 



be expected since the entire Galactic column den- 
sity in this direction is estimatecp] to b e only 



1.9 X 10^° cm-2 (based on |Kalberla et al.||2005 l 



Although both models are statistically ade- 



quate, the Akaike information criterion (Akaike 



1974) suggests one should prefer the model con- 



taining a thermal component. The F test gives a 
null probability of 2 x 10^^^ that the more com- 
plex model would produce such a large improve- 
ment in fit due to chance if the simple power- 
law model were correct. While there are con- 



fusing the HEASARC online calculator, |http : //heasarc 
^nasa . gov/cgi-bin/Tools /w3nh/w3a h~.pl| 
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cerns with using the F test in such a situation 
( [Protassov et al.||2002[ but see also |Stewart]|2009[ ) 
it appears that the statistics do somewhat favor 
a model with a thermal component, although a 
purely non-thermal model cannot be excluded. 

3.2. Orbital variability 

To test for orbital variability, we selected 0.2- 
12 keV source photons from the imaging-mode 
data sets in both our observations and those of 
Homer et al. (2006) and reduced them to the 
solar system barycenter. When combining data 
from multiple instruments, to take into account 
the different orbital coverages and sensitivities, we 
weighted bin values so that all average count rates 
match that observed in the MOSl camera in the 
same observation. For the hardness ratio analy- 
sis, we omitted the PN data so that their different 
energy response and incomplete orbital coverage 
would not skew the results. 

A weighted histogram of these counts is shown 
in Fig. [2] along with a best-fit sinusoid based on 
a finely binned profile. We selected four sinusoids 
as this appears to give a good representation of 
the curve and a time resolution similar to the his- 
togram. A Kuiper test ( Paltani|2004 ) gives a prob- 
ability of 1.3 X 10~^^ (9. Oct) that photons drawn 
from a uniform distribution in orbital phase would 
be this non-uniform; the reduced for a constant 
fit to the histogram is 11.3, with 11 degrees of free- 
dom. While the evidence for variability is strong, 
with only ~ 3 orbits covered by the two observa- 
tions it is not certain that this variability is linked 
to orbital phase, although the fact that the mini- 
mum in the X-ray light curve occurs near orbital 
phase 0.25, when the companion passes closest to 
our line of sight to the pulsar, suggests a link. 
To test this we plotted the photon arrival rates 
for each orbital period separately. Figure [3] shows 
that the flux during each individual orbit appears 
to be lower during phase 0-0.5 than during phase 
0.5-1, which suggests that the variability is indeed 
orbital. Note that Homer et al. (2006) detected 



variability, but having only limited orbital cover- 
age, could not determine whether it was orbital. 
We also computed a hardness ratio (Fig. [2]d), di- 
viding the number of photons harder than 1 keV 
by the number of photons softer than 1 keV and 
comparing the eclipse versus non-eclipse regions 
(for this purpose we defined the "eclipse" region 




1.0 

Orbital phase 



Fig. 2. — Panel (a): Average MOS-equivalent 
count rate in the 0.2-12 keV energy range as 
a function of orbital phase. These data are a 
weighted average of all the imaging-mode data. 
Dashed and dotted vertical lines indicate the be- 
ginning and end of the 1.4 GHz radio eclipse, re- 
spectively. The dashed curve is a best-fit combina- 
tion of four sinusoids. Panel (b): number of MOS 
photons in the 1-12 keV range divided by num- 
ber of MOS photons in the 0.2-1 keV range, also 
as a function of orbital phase. Both panels show 
two cycles for clarity. The companion's closest ap- 
proach to our line of sight to the pulsar happens 
at orbital phase 0.25; that the radio eclipse is off- 
center may be due to eclipsing material being far 
out of the orbital plane, as shown in Figure [5] 
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Fig. 3. — Average count rate in the 0.2-12 keV 
energy range as a function of orbital phase, for 
each orbit independently; the starting MJD for 
each orbit is indicated. The vertical axis is 
MOS-equivalent count rate (based on MOSl/2 
and imaging-mode PN data, scaled to match the 
MOSl count rate), and the dashed horizontal line 
indicates the average count rate over all orbital 
phases and data sets. Where no bar is plotted, no 
data are available. 



to be phases 0-0.5). The reduced for a fit of 
these values to a constant hardness ratio is 8.76 for 
1 degree of freedom, and the probability of such a 
reduced arising if the hardness ratio were con- 
stant is 3.1 X 10"'^ (2.7cr). Thus we see marginally 
significant softening in the eclipse region. 

3.3. Pulsations at the pulsar period 

To test for pulsations, we extracted source and 
background photons from the PN camera using 
the energy range 0.25-2.5 keV, selected to give 
the most significant detection. We barycentered 
the photon arrival times and used the program 
tempcj^and the contemporaneous radio ephemeris 



given in Archibald et al. ( 2009 1 to assign each pho- 



ton a rotational phase. We then tested these pho- 
tons for uniform distribution in phase. The Kuiper 



test (Paltani 2004) gave a (single-trial) null hy- 
pothesis probability of 3.7 x 10~^ (4.5cr) and the 



H test (de Jager et al. 1989) gave a (single-trial) 
null hypothesis probability of 2.4 x 10^'' (4.6(t), 
using an optimal number of sinusoids (two). We 
confirmed that no significant pulsations were de- 
tected with the background photons or with an 
incorrect ephemeris. We also verified that the pe- 
riod predicted by the ephemeris is very close to the 
period at which the significance peaks (holding all 
other ephemeris parameters fixed). 

To estimate the degree to which the X-ray flux 
is modulated at the pulse period, one could simply 
take the lowest bin in the histogram as the back- 
ground level and compute the fraction of photons 
above it; this yields a pulsed fraction of 0.27(9). 
However, this method is subject to large uncer- 
tainties, dependence on binning, and a large sta- 
tistical upward bias due to the fact that we se- 
lected the bin in which signal plus noise is lowest, 
rather than that in which the (unknown) signal is 
lowest. To avoid these problems, we define a root- 
mean-squared pulsed flux by fitting a model F(x) 
with two sinusoids: 



where F is the model mean flux. This is in some 
sense a degree of modulation; if the signal consists 
of a constant background plus a variable emission 



http : //www. atnf ■ cslro . au/research/pulsar/tempo/| 
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Fig. 4. — Panel (a) shows a background- 
subtracted light curve based on 0.25-2.5 keV PN 
photons folded according to the radio ephemeris; 
the sinusoid (drawn with a dotted line) is a two- 
component sinusoid fit to the unbinned photon 
arrival times. Panel (b) shows a 1400 MHz pro- 
file obtained with Arecibo ( Archibald et al.|2009 |. 
Both panels show two cycles for clarity. The 
uncertainty in relative alignment, due primarily 
to the absolute timing uncertainty in the XMM- 
Newton data, is indicated by the horizontal error 
bar. The 20 /is uncertainty relative to the radio 
data due to dispersion measure uncertainty is rel- 
atively unimportant. 



process that drops to zero, this is not directly mea- 
suring the fraction of emission due to the vari- 
able process. (There is a conversion factor that 
depends on the exact shape of the pulse profile.) 
However, this quantity can be computed with sub- 
stantially less uncertainty and bias. Computation- 
ally, we estimate the root-mean-squared amplitude 
in the Fourier domain, based on the amplitudes of 
the two complex Fourier coefficients (since higher- 
order Fourier coefficients are dominated by noise) ; 
since noise always contributes a positive power to 
Fourier coefficients, to reduce the bias we sub- 
tract the expected contribution of noise from the 
squared amplitude of each coefficient before tak- 
ing the square root. We then convert this pulsed 
flux value to a pulsed fraction by dividing by the 
total background-subtracted flux from J1023. 

For the energy range 0.25-2.5 keV we estimate 
a root-mean-squared pulsed fraction of 0.11(2). In 
the two subbands 0.25-0.6 keV and 0.6-2.5 keV, 
we find root mean squared pulsed fractions of 
0.17(5) and 0.14(3), respectively, with profiles that 
are broadly similar and in phase. Above 2.5 keV 
we find no evidence for pulsations, though a Scr 
upper limit on the pulsed fraction is only 0.20. 

4. Discussion 

To summarize our results for J1023, we ob- 
served an X-ray spectrum dominated by a hard 
power-law component, possibly with a small, soft, 
thermal contribution. The emission appears to be 
modulated at the 0.198-day orbital period, with 
substantial dips as the companion passes near the 
line of sight; these dips are accompanied by a pos- 
sible softening in the spectrum. We also found 
that the X-ray emission is very likely modulated 
at the 1.69-ms rotational period of the radio pul- 



4.1. Nature of the X-ray emission 

In a review of X-ray emission from millisec- 



ond pulsars, Zavlin ( 2007 1 describes three primary 



sources of X-ray emission: emission from an in- 
trabinary shock, emission from the neutron star 
itself, which is some combination of thermal and 
magnetospheric emission, and emission from a pul- 
sar wind nebula outside the binary system. As all 
these mechanisms may be operating in J1023, we 
next consider their contributions, if any, to the 
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observed X-rays from this system. 



4.. 1.1. Emission From an Intrabinary Shock 



J1023 likely has a strong pulsar wind, since the 
optical data suggest that the companion is be- 
ing heated by a luminosity of ~ 2Lq from the 
pulsar (Thorstensen & Armstrong 2005). This is 



much greater than the X-ray luminosity of the sys- 
tem (~ O.O3L0), but well below the upper limit 
(~ 8OL0) on the spin-down luminosity of the pul- 
sar ( Archibald et al. 2009 1 . If material is leaving 



the companion, either through Roche-lobe over- 
flow or a stellar wind, we should expect an in- 
trabinary shock, where this flow meets the pulsar 
wind. Such a shock could readily produce power- 
law X-ray emission ( |Arons fc Tava^[l993| . If lo- 
calized, it could easily account for the orbital mod- 
ulation we observe in the X-ray emission. If the 
emission were due to Roche-lobe-overflowing ma- 
terial meeting the pulsar wind we might expect it 
to be localized at or near LI. Given the system 



geometry described in Archibald et al. (2009) (a 
46° inclination angle, corresponding to a pulsar 
mass of I.AMq, and a Roche- lobe-filling compan- 
ion; see Figure [5|, the LI point itself is eclipsed 
by the companion for 0.32 of the orbit, centered 
on orbital phase 0.25. We would expect spectral 
softening during this period, since the relatively 
hard power-law emission is blocked but the, pre- 
sumably softer, emission from the neutron star is 
not. A larger X-ray-emitting region, due either to 
material streaming away from the LI point or to 
a stellar wind shock, would result in a broader, 
shallower, and potentially asymmetric dip in the 
X-ray light curve. 

The possibility of a wind from the companion 
is interesting, since it is not clear how Roche-lobe 
overflow would provide the ionized material caus- 
ing the radio eclipses, which occur well above the 
orbital plane, roughly centered on the companion's 
closest approach to our line of sight. Archibald] 
et al. (2009) reported dispersion measure (DM) 
changes of ^ 0.15 pc cm"-^ above a cutoff fre- 
quency of ^ 1 GHz at radio eclipse ingress and 
egress. Assuming that this is the plasma frequency 
implies an electron density of ~ 10^° cm~^; com- 
bined with the excess DM measurement, this im- 
plies a layer of thickness ~ 4 x 10^ cm, rela- 
tively thin compared to the orbital separation of 
1.2 X 10^^ cm. This thin sheet of material could 
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Fig. 5. — J1023 system geometry at phase 0.25 as 
seen from within the plane of the sky. The system 
orbit is in and out of the page, the dashed line 
is our line of sight to the pulsar, and the dotted 
line is our line of sight to LI, which is eclipsed 
by the Roche-lobe-filling companion. The pulsar 
is indicated by a not-to-scale circle, but the com- 
panion is drawn to scale. The tick mark on the 
line connecting pulsar and companion marks the 
center of mass of the system. This system geom- 
etry assumes a pulsar mass of 1.4 Mq; a higher 
mass would mean a larger system seen more nearly 
face-on, but the range of plausible angles is con- 
strained to ^ 53°-34° by the narrow range of plau- 
sible pulsar masses 1.O-3.OM0 (Archibald et al 



2009 ) . Note that at this orbital phase the X-rays 
are near their minimum, and low radio frequen- 
cies from the pulsar are eclipsed, indicating the 
presence of ionized material well out of the orbital 
plane. 
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in principle arise as the shock front where the pul- 
sar wind meets a wind from the companion, but 
it is difficult to explain the companion's wind be- 
ing strong enough to support such a shock: given 
the line of sight geometry of the system, the ma- 
terial causing the radio eclipses must be about as 
far from the companion as it is from the pulsar. 
If we assume that the highly relativistic pulsar 
wind provides the ~ required by the compan- 



ion heating models of [Thorstensen fc Armstrong] 
( 2005 1 , it seems difficult to explain how the com- 



panion could have a wind of comparable pressure 
anywhere along the line of sight. Thus the radio 
eclipses remain a mystery, and Roche-lobe over- 
flow seems a more likely explanation for the origin 
of the intrabinary shock. 

4- 1.2. Emission From the Neutron Star 

Some millisecond pulsars, binary or isolated, 
have X-ray emission modulated at the rotational 
period. This emission is some combination of mag- 
netospheric emission, from high-energy particles 
in the magnetosphere, and thermal emission, from 
polar caps heated by bombardment by high-energy 
particles moving along the open field lines. Ob- 



servationally, Zavlin (20071 draws the distinction 



that magnetospherically dominated emission has 
high pulsed fractions, narrow peaks, and power- 
law spectra, while polar-cap-dominated emission 
has lower pulsed fractions, broader peaks, and 
thermal spectra. 

tend 



Zavlin| also suggests that pulsars 
with spin-down luminosities > 10'^^ 



erg s 



to be magnetospherically dominated while those 
with spin-down luminosities < 10'^^ erg s^^ tend 
to be thermally dominated, but since few examples 
are known, this classification is tentative, and the 
spin-down luminosity from J 1023 is only known 



to be < 3 X 10^^ erg s"^ (Archibald et al. 20091 



In any case, classification of the pulsations we ob- 
serve in J1023 is difficult. 

In systems in which thermal emission from the 
neutron star provides all the detectable X-rays, 
detailed modelling indicates that the pulsed frac- 
tion should normally be < 50% (that is, pulsed 
emission must be accompanied by roughly equal 
or greater unpulsed emission) due to light bend- 
ing and the large size of the polar caps, though for 
certain combinations of parameters higher pulsed 
fractions can occur ( Bogdanov et al.|["2008 ). If the 
pulsations are purely thermal, then our spectral 



upper limit on the thermal fraction of 0.06(2) is 
difficult to reconcile with the fact that we observe 
a fractional modulation of 0.11(2) at the pulsar 
period. Thus it seems likely that at least some 
magnetospheric emission is present, as it has a 
nonthermal spectrum and can more readily have 
higher pulsed fractions. The poor signal-to-noise 
in our observations, due to the high background 
(from the system itself as well as instrumental) 
and the scarce photons, makes it impractical to 
determine the sharpness of the pulse profile or the 
hardness of the pulsations. In any case the pul- 
sations, if real, are a clear sign of X-ray emission 
from the pulsar itself. 

4- 1.3. Emission From a Pulsar Wind Nebula 

Nebular emission from pulsars arises when the 
particle wind driven by the pulsar's magneto- 
spheric activity flows out of the pulsar's immediate 
neighborhood and meets the surrounding medium; 



for reviews see Gaensler & Slane (2006) or Kaspi 



et al. ( 2006 ) . Pulsars like J1023 whose wind is con- 



fined by the ambient interstellar medium exhibit 
nebular emission that generally takes a cometary 
form, with an arc-like bow shock preceding the 
pulsar and a "trail" of ejected material streaming 
back along the pulsar's track. 

The angular size of bow shock emission depends 
on the pulsar's spin-down luminosity, its proper 
motion, and the local interstellar medium den- 



sity. Kargaltsev & Pavlov (20081 give a formula 



for predicting the "stand-off" angle of the X-ray 
bow shock from the pulsar: 



= 5.4" 



,-1/2,-1^.1/2^-2 



' — i pi 

'0.1 f^l9 ^35 -^1.3' 



where no.i is the mean density of the interstellar 
medium divided by 0.1 cm~^, ^ig is the proper 
motion divided by 19 mas yr~i, £^35 is the spin- 
down luminosity divided by lO'^^ erg s~^, and 
-D1.3 is the distance to J1023 divided by 1.3 kpc. 
We have supplied best-guess parameters for J1023 
( Archibald et al.|[2009 ), so that absent further in- 
formation, we might expect a stand-off distance of 
^ 5" for the X-ray bow shock. Thus although neb- 
ular emission is not resolved in our XMM-Newton 
images, higher-resolution images with the Chan- 
dra X-ray Observatory, or in the Ha or radio 
bands, might yet resolve a bow shock. 



Kargaltsev & Pavlov (2008) find that the ra- 



tio of neutron-star emission to nebular emission 
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generally lies between ~ 0.1 and 10, which sug- 
gests that nebular emission may produce some 
of the X-rays we observe. Moreover, |Kargaltsev| 
Ifc Pavlovl find that nebular emission from bow 
shocks has a power-law spectrum with photon in- 
dices 1 < r < 2 ( [Kargaltsev fc Pavlov||2008[ ), con- 
sistent with what we observe from J1023, though 
the latter is somewhat harder than most pulsar 
wind nebulae. The orbital variability, on the other 
hand, argues against that component of the emis- 
sion arising from an extrabinary pulsar wind neb- 
ula. 

4.2. Comparison to similar systems 

In quiescence, most LMXBs have spectra that 
are dominated by thermal emission from the neu- 
tron star. This emission is typically consistent 
with a hydrogen atmosphere model with effective 
radius 10 km (iBrown et al.|1998|). Some LMXBs 



in quiescence also have a small power-law compo- 
nent of photon index 1 < F < 2 in their spec- 



trum (Campana et al. 1998), while a few quies- 



cent LMXBs have spectra completely dominated 
by this non-thermal component. For example, the 
quiescent LMXB SAX J1808.4-3658 has a spec- 
trum that is a fairly hard power law with no de- 



tectable thermal emission (Heinke et al. 20091. 
J 1023 resembles this latter category, although if a 
thermal component is present its effective radius 
appears to be substantially less than 10 km and 
the power-law spectral index is fairly hard. On 
the other hand, J1023's radio and possible X-ray 
pulsations distinguish it from all known quiescent 
LMXBs. 

SAX J1808.4-3658 is of particular interest be- 
cause it is the prototype and best-studied exam- 
ple of a class of LMXBs from which millisecond 
X-ray pulsations have been detected during ac- 
tive phases. Many authors have suggested that 
SAX J1808.4— 3658 should turn on as a radio pul- 



sar during quiescence (e.g. 


Wijnands & van der 


Klis 19981 


Chakrabarty & Morgan 


1998 




Cam- 


pana et al. 


2004), but no radio pulsations 


have 



been detected in spite of thorough searches (e.g. 
Burgay et al]|2003| [Tacolina et al]|2009[). The 



non-detection of radio pulsations could be be- 
cause SAX J1808.4— 3658 is shrouded by previ- 
ously ejected ionized material, or it could sim- 
ply be an issue of unfortunate beaming. In any 
case, it is natural to compare SAX J1808.4— 3658 



to J1023. 

During its active phases, SAX J1808. 4-3658 
reaches a sustained luminosity of ^ 4x 10^^ 



erg s 



(Galloway & Gumming 2006), in sharp contrast 
to the upper limit of 2 x lO"^"* erg s" 



available 



for J1023 during its disk phase (Archibald et al 



2009) 



-3658 



In its quiescent state, SAX J1808.4 
shows X-ray emission dominated by a hard power 
law (F = 1.74(11)), with a small (< 0.1) frac- 
tion of thermal emission, and a t otal X-ray lu- 
minosity of 7.9(7) x 10'^^ erg (Heinke et al. 
|2009J . This is somewhat similar to what we ob- 
serve in J 1023, which has a total luminosity of 
9.9(5) X 10'^^ erg s~^, although the spectral index 
of 1.26(4) is somewhat harder than that seen from 
SAX J1808.4-3658. If J1023 has a thermal com- 
ponent, its power-law is substantially harder, with 
a spectral index of 0.99(11). The thermal compo- 
nent would have an effective radius of 0.7^q'^ km 
and supply 0.06(2) of the total luminosity. While 
the spectral index seen in J1023 is harder than 
that seen from SAX J1808.4-3658, it is in fine 
with the hard spectral indices seen from other 
X-ray-detected radio millisecond pulsars (for ex- 
ample 47 Tuc W, though PSR J1740-530 is softer; 
see below). Since SAX J1808.4-3658 is known 
to cool quickly, even thermal X-ray pulsations in 
quiescence would be evidence of magnetospheric 
activity heating its polar caps. [Gampana et al.| 
(2002) searched for such X-ray pulsations from 



SAX J1808.4— 3658 in quiescence, but were unable 
to conduct meaningful searches due to its faint- 
ness; SAX J1808.4— 3658 is roughly 3.5 kpc away 



(Galloway & Gumming 2006) compared to the 



roughly 1.3 kpc we assume for J1023 (Archibald 
etar][2009l ). 



The similarity of the X-ray properties of SAX 
J1808.4— 3658 in quiescence to those of J1023 sug- 
gests that SAX J1808.4-3658 may indeed har- 
bor a radio millisecond pulsar in quiescence; con- 
versely, the same similarity suggests that J1023 
may undergo episodes of active accretion. 

Among radio pulsars, the two best-studied ex- 
amples which closely resemble J1023 are 47 Tuc W 
( IBogdanov eVal] |2005| and PSR J1740-5430 
( D'Amico et al.|2001 ). Both pulsars exhibit large, 
variable radio-frequency eclipses, and both are 
thought to have somewhat massive unevolved 
companions (> O.13M0 and 0.19-0.8 Mq, re- 
spectively). For comparison, J1023's companion 
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is thought to be ~ 0.2 Mq, and shows a spectrum 



characteristic of a main-sequence star ( 


Archibald 


et al.|2009 


Thorstensen & Armstrong|2005 


I. The 



radio eclipses are in all three cases evidence for 
material leaving the companion and presumably 
being expelled from the system, and in all three 
cases the companion appears to fill, and perhaps 
overflow, its Roche lobe (Bogdanov et al. 2005 
Ferraro et al.| |2001| [Thorstensen fc Armstrong] 
20051). 



Bogdanov et al. (2005) studied 47 Tuc W and 



found an X-ray spectrum consisting of a domi- 
nant hard power-law (F = 1.14(35), contributing 
~ 75% of the total flux) plus a thermal compo- 
nent. They also observed orbital variability, and 
in particular a substantial X-ray eclipse spanning 
phase^ 0.15-0.45 and centered slightly ahead of 
the optical minimum, at phase 0.2. This is ac- 
companied by substantial softening of the X-ray 
spectrum during eclipse. Given these observa- 
tions, ,Bogdanov et al.| suggest that the power- 



law X-ray emission originates from a shocked re- 
gion where material overflowing from the compan- 
ion of 47 Tuc W meets the pulsar wind and is 



blown out of the system. Bogdanov et al. (2010) 
studied Chandra X-ray Observatory observations 
of PSR J1740-5430. Ahhough limited by the 
scarcity of photons, they found a spectrum con- 
sistent with a power law of index F = 1.73(8), or 
with a somewhat harder power-law plus a black 
body component. Their estimates of variabil- 
ity were limited by restricted orbital coverage as 
well as by a paucity of photons, but they found 
marginal evidence for orbital modulation, in par- 
ticular a possible decrease in luminosity during the 
radio eclipses, possibly accompanied by a soften- 
ing. Our observations of J1023 are more photon- 



starved than those of Bogdanov et al. (2005) of 



47 Tuc W, but we do see evidence for eclipses, 
and possibly for softening during eclipses. In the 
case of J 1023, the possibility of overflowing gas is 
supported by the recent disk phase. 

The similarity of the observational properties 
of these three millisecond pulsars suggests that 
their companions are all overflowing their Roche 
lobes, so that mass flows through the LI point. 



'''Note that these authors use a diflferent convention for or- 
bital phase, setting phase 0.5 to the optical minimum, while 
we set phase 0.25 to the optical minimum. All phases given 
here have been converted to our convention. 



where it encounters the pulsar wind and is swept 
away. On the other hand, the similarity also sug- 
gests the possibflity that like J1023, 47 Tuc W or 
PSR J1740— 5430 may undergo episodes in which 
the mass transfer rate from the companion in- 
creases and an accretion disc forms. Such episodes 
would presumably be signalled by optical changes, 
possibly by extinction of the radio pulsations, and 
presumably by modest X-ray brightening. All to- 
gether, these observations suggest that as a system 
reaches the end of its life as an LMXB, its accre- 
tion rate drops and the millisecond pulsar becomes 
active. After this point, the wind from the pul- 
sar generally prevents mass leaving the compan- 
ion from forming an accretion disk or entering the 
pulsar magnetosphere. Temporary increases in the 
mass accretion rate may allow the occasional sup- 
pression of the wind and formation of an accretion 
disk, as was observed in J1023 in 2001. 

5. Conclusions 

J1023 shows X-ray emission probably consisting 
primarily of emission from an intrabinary shock, 
plus a smaller amount of emission from the pul- 
sar itself and possibly some unresolved nebular 
emission. The shock emission shows substan- 
tial variability that appears to be linked to or- 
bital phase, as well as possible spectral soften- 
ing during eclipses. Similarities to 47 Tuc W and 
PSR J1740-5340 suggest that this may be due 
to gas overflowing from the companion meeting 
the pulsar wind in a shock close to the compan- 
ion. The emission from the pulsar itself, if real, 
is probably due to some combination of curva- 
ture radiation in the magnetosphere and thermal 
radiation from polar caps heated by high-energy 
particles streaming downward, though which ef- 
fect dominates is unclear. Further observations to 
provide improved spectra should help clarify the 
origin of both pulsed and unpulsed emission, as 
would further high-time-resolution observations to 
allow better measurement of pulse profiles, pulsed 
hardness ratios, and orbital variations of pulsed 
fraction. X-ray observations with higher spatial 
resolution might also resolve extended nebular 
emission around J1023. Better understanding of 
this object promises to clarify the transition from 
LMXB to radio millisecond pulsar, and may show 
that other known systems are currently in this fas- 
cinating transition state. 
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